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Abstract 

It is proposed that soft particle production in hadronic colhsions 
is dominated by multiple gluon exchanges between partons from the 
colliding hadrons, followed by radiation of hadronic clusters from the 
coloured partons distributed uniformly in rapidity. This explains nat- 
urally two dominant features of the data: (a) The linear increase of 
rapidity spectra in the regions of limiting fragmentation and, (b) the 
proportionality between the increasing width of the limiting fragmen- 
tation region and the height of the central plateau. 
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1. Recently, a substantial evidence is accumulating that particle produc- 
tion in hadron-hadron, hadron-nucleus and nucleus-nucleus collisions at high 
energy satisfies the principle of limiting fragmentation H!] in a much wider 
range of rapidities than originally proposed. This effect, first seen in pp j2] 
and in p- Emulsion collisions [3] , was recently studied in ci — Au and Au — Au 
interactions by the PHOBOS collaboration jHEl- 




Figure 1: Particle density in pseudo-rapidity plotted versus Y — rj. (a) p — p 
collisions j2]. (b) nucleus- nucleus collisions 0. Lines are drawn to guide the 
eye. 

The phenomenon is illustrated in Fig. 1, taken from [2] (la) and j3] 
(lb), where particle density in pseudo-rapidity is plotted versus the difference 
between the beam rapidity Y and pseudorapidity rj of the particle. The data 
of UA5 collaboration |2] on pp collisions and recent data from PHOBOS 
collaboration j3] on Au-Au collisions are shown. One sees three prominent 
features, common for the two data sets: (i) Except at very small Y — rj, 
particle density in the fragmentation region increases linearly with increasing 
Y—rj] (ii) This linear increase is followed by a "plateau" in the central rapidity 
region; (iii) The width of the central plateau grows with increasing Y only 
very slowly, if at all. This last feature implies that with inreasing Y (i.e. 
increasing energy) the range of the limiting fragmentation region increases 
proportionally to Y . 
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It should be emphasized that these features are in blatant disagreement 
with the principle of boost-invariance in particle production [Uj. They are 
thus difficult to understand in the standard description of production pro- 
cesses where one expects the particle density to be dominated by the central 
plateau (HI E], as is common in many current models. In contrast, as is seen 
in the Fig. 1, the "central plateau" occupies only a fraction of the available 
rapidity range. 

In the present note we show that all these features can be understood in a 
picture where particle production proceeds by a number of colour exchanges 
between the two sets of partons one from the projectile and one from the 
target. These colour exchanges lead to creation of the colour charges which 
emit the observed particles by the bremsstrahlung process. If the original 
partons in each of the colliding hadrons are uniformly distributed in rapidity 
(i.e. if they satisfy the Feynman dx/x rule) the resulting distribution of 
observed particles is linear m. Y — y. Noting that only partons with life- 
time longer than the time tq, needed for the colour exchange to take place, 
can participate in the process, we conclude that this linear increase of the 
spectrum must stop at a rapidity y = yo, depending on the parton transverse 
mass /i and on tq. Thus the linear increase is followed by a plateau for y 
smaller than y^. By postulating that this picture is valid in the cm. frame 
of the collision, thus violating explicitely the boost-invariance, one accounts 
for the gross features of the data. 

In Sections 2 and 3 a more detailed description of the model is presented. 
In Sections 4 and 5 the consequences of the model for particle spectrum are 
described. Discussion and conclusions are given in the last section. 

2. We follow the standard approach to multiparticle production thus 
accepting that high-energy collision of two hadrons can be described by colour 
exchange between the partons from the projectile and from the target jH]^. 

The new idea which we propose in this note is that this mechanism is 
realized by colour exchange between several pairs of partons chosen at ran- 
dom, one from the projectile and one from the target. The members of each 
pair moving in opposite direction, each one radiates the observed particles 
(or particle clusters^) in the process of bremsstrahlung ^Tj. Note that by 
this postulate we abandon the Feynman assumption ^ that the interaction 
is dominated by "wee" partons, i.e. by partons with small rapidities. 

^There are many models of this type and many reviews. See, e.g. pH llOj . 
^E.g.in the form of hadronic resonances. 
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To justify this idea we observe that it can be reahzed simply by one gluon 
exchange process. Indeed, the gluon being particle of spin one, its exchange 
gives the energy-independent cross-section. This means that the probability 
of interaction between two partons by exchanging a gluon is independent of 
their relative rapidity. In short, our assumption is justified by the existence 
of vector particles in QCD. 

To obtain specific predictions one needs to know the shape of the parton 
distribution. We simply accept the standard idea that the partons in each of 
the colliding hadrons are distributed uniformly in rapidity and we can thus 
write their distribution in the form [TT] 



where b is the parton density per unit of rapidity and z+ = {E + pL)/{Ei + 
Pi), with {E,pl) are the energy and longitudinal momentum of the parton, 
whereas {Ei, Pi) are the energy and momentum of the beam^. 

This picture makes sense only if the partons in one projectile can be 
treated as independent from those in the other one. This may be justified 
if the rapidity separation is large enough. When the rapidity separation is 
small, however, (i.e. in the region close to the rapidity of the center of mass) 
the partons from the two projectiles mix up and one cannot expect them 
to act independently. Also the notion of the colour separation looses the 
meaning. This case thus demands a special attention. 

This "very central" region is best studied in the overall cm. frame. In this 
frame the parton energies are not large and thus their distribution rapidly 
fluctuates in time. To participate in the collision process, however, the life- 
time of a parton must be substantially longer than the time tq needed for the 
interaction to take place. This condition allows to estimate the effective size 
of the rapidity region which does not contribute to the particle production. 

To see this we observe that the life-time of a high-energy parton with 
transverse mass equal to fi can be estimated from the uncertainty principle 
as 



■^These formulae apply for the right-moving system. For the left-movers one should 
replace z+ by z- = {E -pl)/{Ei - Pi). 




(1) 



r 



7//i = E/^^^ 



(2) 
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where 7 is the Lorentz factor. From the condition r ^ tq we obtain /ie^'/2yU^ ^ 
To where y' is the rapidity of the parton. This imphes 

2 

ey' > » 2/xro; ^ > » (3) 

One sees that the condition ^ restricts substantially the rapidity of 
partons which can participate in particle production, as one may expect tq 
to be of the order of 1 fermi(ro ~ 1/pt, where pt is the transverse momentum 
exchanged in the interaction)^. 

3. The emission of particle clusters in the bremsstrahlung process was 
analyzed some time ago by Stodolsky jjjj . We follow his approach and write 
the particle distribution in the form 

dN{x+) = a{l - x+y-^^ (4) 

where a is the density of emitted hadrons per unit of rapidity and x± = 
(e ± ql) I [Ei + Pj), with (e, g^) being the energy and longitudinal momentum 
of the emitted cluster. 

Denoting by A the fraction of "active" partons, i.e., the partons which par- 
ticipated in the colhsion and using ((H), the distribution of the bremsstrahlung 
products is 



J z Zx- 



a 1 



X-\- \ (lX-\ 



X4 



(5) 



where the lower limit of integration z is 

z = 'max{x+, Zq) (6) 

with Zq determined by the energy below which a parton does not live long 
enough to undergo a soft interaction and therefore also does not radiate(c.f. 
and the related discussion in the previous section). 
By changing the variables: 

u = 1 — x+; u = 1 — z; z^ = 1 — ut (7) 



^This restriction is of course much less effective for hard collisions where the interaction 
time may be very short. 
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we obtain 

rIN rIN r^^/u 

oiN^ = ^ = Aa6M"+^-^ / dtil - ut)-%l - tY~H'-' (8) 
dx+ ay Jo 

4. We are mostly interested in the particle distribution for rapidities out- 
side the projectile fragmentation region, i.e, for small x+ ~ 0. The formula 
(jHI) shows that we have to consider two cases. 

For x+ < zq we have u = 1 — zq. Thus, in the limit x+ — we obtain 

— ^ Xab 1'^-^ (9) 

dy Jo 1-t ^ ^ 

One sees that the result is independent of m = 1 — a;+, i.e. we obtain a plateau 
for y < yo. 

For X > ZQwe have i = x+, i.e. u = u. Consequently, (jHl) can be rewritten 

as 



dy r(a + b) 

where F is the hypergeometric function. Note that the result is perfectly 
symmetric with respect to a and h. 

To see the behaviour at a;+ ~ 0, i.e. u ^ 1 we use the formula giving 
expansion of -F(a, 6; a + 6; u) around u = 1 ^21- In the limit of small Xj^ this 
gives 

^ = \ah[2^{] + l)-^{a + ])-^{h + 3)-\ogx+] = 

= \ab [2^(1) - ^(a) - ^(6) + log(M/m) +Y -y] (11) 

where m is the transverse mass of the emitted cluster and M is the mass of 
the incident particle. We thus obtain a linear increase with increasing Y — y, 
as observed in the data. 

The behaviour in the fragmentation region x+ ~ 1, u ~ 0, is best seen 
from(|TU|). The result is 

*v ^ /(« + i)r(^ +i) (,, _ (12) 

dy T[a + b) 

which shows a deviation from the linear increase unless a + b = 2. 
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5. The distribution discussed in Section 4 can be justified for positive 
rapidities but is not applicable in the negative rapidity region. This is the 
standard problem in the bremsstrahlung model JT]. The reason is clear: 
Eq (jU was obtained by requiring conservation of the sum (e + ql), ignoring 
entirely conservation of the difference (e — ql) (i.e. ignoring entirely the 
target). 

Since the division into positive and negative rapidities is frame-dependent, 
the region of applicability of (jH)) is also frame-dependent. To fix this, we are 
again forced to use the hypothesis that our considerations are vahd in the 
cm. frame of the collision. Of course any other frame boosted by less than 
Ho is equally good. 

In the actual calculations (shown in Fig. 2) we cut the distribution for 
negative cm. rapidities using a simple prescription to multiply the distribu- 
tion (jH)) by the correcting factor 

x+ + X- e^y + 1 

One sees that for positive (large) y the correction factor is unimportant. 
On the other hand, it cuts exponentially the distribution for negative y. A 
similar procedure must be, of course, applied also to the other projectile, 
where one takes $_(?/) = $+(—?/) = 1 — $+(?/)• The observed distribution is 
the sum of two contributions, one from the right-moving system and another 
one from the system moving to the left. 

Since the model predicts a plateau in rapidity between — ?/o and +?/o, the 
exact form of the cut-off is not essential, as long as it is ineffective beyond 
the plateau region [TT|. 

In Fig. 2 the calculated dN/dy is plotted versus Y — y ioi a = 1, Xb = 1, 
yo = 2, M = m, two energies and three values of the parameter b. One sees 
that the numerical results confirm the semi-quantitative conclusions given in 
the previous section. One also sees that they resemble nicely the data shown 
in Fig. 1. 

6. Several comments are in order. 

(i) The idea that the colour charges created in the first step of the col- 
lision are responsible for particle production is rather general and may be 
implemented in many ways. An interesting possibility is to consider a more 
detailed model assuming that the active (radiating) partons are colour octets 
and the density of hadronic clusters is proportional to the local density of the 
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Figure 2: Particle density calculated from (jH)). Parameters as shown in the 
figure. The dashed lines show the effect of the cut-off (fT3|) 

chromoelectric field created by these partons. For a given rapidity y there 
are, say, octets moving to the left (i.e. having rapidities smaller than y) 
and colour octets moving to the right (i.e. those with rapidities greater 
than y). Colour conservation implies that the representation formed by 
left movers is conjugate to the representation Rr formed by right movers: 
Rr = Rl 

If the local energy density S{y) of the chromoelectric field is proportional 
to the quadratic Casimir operator C2 (Rl) = C2 (Rr) and for < tir 

i?^ = 8 (g) 8 (g) . . . 8 (n^times) (14) 

then the average energy density 

(^)ni,n. = ^oCAmin {ul, Ur) (15) 

where Ca = 3 is the value of the Casimir operator for the adjoint represen- 
tation and Sq is a constant. The linear increase of hadron rapidity density is 
obtained if the rapidity distribution of radiating partons is uniform in both 
hemispheres. The maximum of the single particle distribution is obtained for 
the rapidity corresponding to rii = ur. In p — p collisions this is of course 
the rapidity of the center of mass. 

(ii) We have worked out in detail the hypothesis that the production 
of final hadrons from the colour charges proceeds by the bremsstrahlung 
process. This formulation is by no means unique. Conclusions similar to 
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the ones reached in this paper can be obtained if the production is described 
by breaking of colour strings spanned between a parton from the projectile 
and a parton from the target. In this case the cut-off (|T^ is not needed. 
Instead, to obtain the plateau in the central region one may postulate that 
the difference between the ends of a string contributing to particle production 
must exceed 2?/o- This does not change the main results of the model. The 
essential point, needed to obtain the linear increase of the rapidity spectrum, 
is the flat distribution (in rapidity) of the radiating partons and the fiat 
distribution of clusters in string decay (with energy independent density). 

One should note, however, that these two versions of the model give ob- 
servable differences for asymmetric heavy ion collisions (in particular p-A and 
d-A collisions) and for forward-backward correlations in particle multiplici- 
ties. It may thus be interesting to study these correlations experimentally. 

(iii) The model contains several unknown parameters which, fortunately, 
have a well-defined physical meaning. Their determination from the data 
may thus give an interesting insight into structure of hadrons relevant for soft 
interactions. For example, a determination of the effective parton density at 
low momentum transfers (described by the parameters b and A) is of clear 
interest. When applied to nuclear collisions, this would allow to investigate 
the relation between the effective parton densities and other parameters such 
as the number of wounded nucleons JSj and/or number of collisions. 

(iv) Our argument explains only the gross features of the data. To obtain 
a more detailed description, it is necessary to include at least the effects of 
cluster decays^. This seems feasible, particularly in pseudorapidity, where 
the isotropic clusters have a well-known decay distribution. 

7. In conclusion, to explain the observed strong violation of boost in- 
variance in rapidity spectra, we have proposed a two-step mechanism for soft 
particle production in hadronic collisions. The first step is the multiple gluon 
exchange between the partons from the two colliding hadrons. In the second 
step, partons which were involved in this process radiate hadronic clusters. 
This mechanism provides a natural explanation of the observed rapidity spec- 
tra, in particular their linear increase with increasing rapidity distance from 
the maximal rapidity. Also the short plateau in the central rapidity region is 
naturally obtained. The scheme can be applied also to nuclear collisions and 
is flexible enough to account for the gross features of the data. All parame- 

^Also the change of variables {y <-> 7j) may be important, particularly in the region 
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ters needed in this description have a well-defined physical meaning and thus 
their determination from the data would give useful information on hadron 
structure in the non-perturbative region. 
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